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ABSTRACT 


The effect of phosphorus as a modifier on titania supported chromium oxide was studied for 
the oxidative dehydrogenation (ODH) of propane. A series of catalysts of unmodified and 
phosphorus modified chiomia/titania in three different chromium to phosphorus molar 
ratios of 2:1, 1.5:1 and 1:1 were prepared by the incipient wetness impregnation method. 
The unmodified and modified chromium oxide catalysts were characterized for their surface 
area and by XRD, EPR and TPR techniques. No crystalline peak of or related to Cr 203 , and 
P2O5 were detected in the XRD of prepared samples. Clustered Cr^’^-O-Gr^^ species in 
addition to the y-Cr^^were detected by EPR. From the TPR studies, it was found that the 
reducibility (Tmax) of the surface chromium species was not affected for the modified 
catalysts. The H/Cr molar ratio of the samples, however, decreased with phosphorus 
addition. The ODH of propane reaction studies revealed the effect of phosphorus on the 
propene selectivity and propene yield of the surface chromium oxide species. The 
conversion progressively decreased with phosphorus addition. Propene selectivity at iso- 
conversion increased till a chromium to phosphorus ratio of 1.5 and then decreased. Two 
power law (PL) models were considered to explain the reaction data. The power law model 
that considers a mechanism in which the propene is a primary product and the carbon 
oxides are secondaiy products is observed to be the best model to fit the reaction data. The 
kinetic parameters are obtained by minimization of an objective function using Genetic 
Algorithm. The pre-exponential factors decreased with phosphorus addition for most of the 
samples, however, the relative decrease was different. The activation energies and partial 
pressure exponents are not affected or affected randomly. Based on estimated kinetic- 
parameters the propene yield was calculated and observed to be maximum for the 
1.5CrlPTi catalyst. For this catalyst the ratio of the rate constant for propene degradation 
(k 2 +ki) to propane ODH (ki) was the minimum. 
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Chapter-1 


Introduction 

Light alkenes are important building blocks in the modem petrochemical industry. They 
find application for the production of a wide variety of products. Consequently, the 
demand for alkene has been increasing. Due to this rising demand of alkenes, conversion 
of alkanes into their corresponding alkenes is an industrially significant process. Low 
cost of light alkanes and the fact that they are generally environmentally acceptable due 
to their low chemical reactivity have made them the ideal feedstock for alkene 
production. Accelerating global demand for alkene has spurred substantial interest in the 
development of alternative routes other than the traditional processes. Presently, steam 
cracking, catalytic cracking and dehydrogenation are mainly used as the direct route for 
the conversion of lighter alkanes to high value olefins. 

The propene demand, mainly for production of polypropene, acrolein and acrylic 
acid, is increasing at a higher rate than the demand for ethylene and other light olefins co- 
produced in the traditional processes such as steam cracking and FCC of natural gas and 
oil fractions. Since the last few years, catalytic dehydrogenation of propane is also a topic 
of commercial as well as scientific interest as the better understanding of these reactions 
allows exploitation of cheap and abundant resources and selective production of propene. 
Despite the apparent simplicity of the reaction, however, the dehydrogenation processes 
possess several thermodynamic constraints, such as high endothermicity which require 
severe conditions of high temperature and low pressure [1]. Thermal cracking side 
reactions favored by high temperatures, significantly affect the stability of propylene and 
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alkenes produced. Thus, coke is deposited on all dehydrogenation catalysts, causing a 
drop in their activity [2, 3], Consequently, frequent regeneration and complicated reactor 
designs are required. The thermodynamic and kinetic constraints have stimulated the 
search for better solutions. 

To meet with all the above-sited challenges a potential new route, Oxidative 
Dehydrogenation (ODH), has been introduced for the production of alkenes [1, 4]. With 
the successful development of ODH, high olefin yields will be possible through the 
conversion of much smaller volumes of alkanes. Oxidative dehydrogenation (ODH) of 
alkanes is thermodynamically favored over a wide temperature range and offers the 
potential for the synthesis of alkenes with significant savings in energy and feedstock 
costs [1, 4-7]. 

Oxidative dehydrogenation (ODH) overcomes the thermodynamic limitations 
faced in the direct catalytic dehydrogenation process and the rapid coking of these 
dehydrogenation catalysts resulting in short catalyst life. In the ODH reaction, hydrogen 
is oxidized to water and there is no equilibrium limitation. 

CsHg + 02 = CsHe .+ H 2 O AH=-86kJ/mol (1.1) 

This process has the potential to overcome the major technical problems associated with 
catalytic dehydrogenation. However, other problems arise that are related to the removal 
of the heat of reaction, control of selectivity due to the formation of undesired by- 
products and carbon oxides, the flammability of reaction mixture and the possibility of 
the run-away of the reaction. Selectivity to olefins remains a serious problem, as it limits 
the maximum achievable yield. Therefore, developing means of maximizing the alkene 
formation by minimizing side reactions is the major issue for the oxidative 
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dehydrogenation reactions. This can be achieved by the design of a catalyst that 
kinetically controls the oxidative dehydrogenation reactions to optimize the selectivity at 
higher conversions. 

Recent studies have focused on using supported metal oxides for the ODH 
reactions due to their high degree of mechanical strength, better thermal stability and 
larger surface area [8, 9]. Previous studies dealing with supported metal oxide catalysts 
have shown that the nature of the surface metal oxide species depends on the oxide 
support, metal oxide loading, preparation method, and the presence of additives [10], Of 
the various supported metal oxide systems supported chromium oxides are effective 
catalyst for ODH of propane [11-18]. They have been successfully used for 
dehydrogenation and other industrially important reactions [19-21]. Catalytic properties 
of the surface chromium oxide (chromia) species are observed to change due to the 
interaction between support and the chromium oxide species. It is also observed that the 
molecularly dispersed chromia species are more active for ODH of propane compared to 
bulk Cr 203 [22]. It has been found that the specific activity of the surface chromium 
oxide species is highest for Cr 203 /Ti 02 compared to chromium oxide supported on other 
different supports (Si02, Si02-Al203 and AI2O3) [23]. However for Cr203/Ti02 the 
propene selectivity was small. Therefore, it is essential to design a catalyst that optimizes 
both conversion and selectivity of ODH reaction. 

Selective oxidation of hydrocarbons, of which ODH is a special case, requires 
redox and Lewis type acidic properties [24]. Furthermore, an acid-base bifunctional 
catalyst possessing weak acid- base pairs exhibits better catalyst activity and selectivity 
for the ODH reactions [25]. Thus, an ideal catalyst would require a proper combination of 
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redox and acid-base properties. Changes in the redox and acid-base properties can be 
brought about by using additives/modifiers. Such studies over Cr203/Ti02 catalyst were 
performed by Cherian et al, [ 26 ], Two types of additives were observed-interacting 
additives (Na) that co-ordinate with the surface chromia site and non-interacting additives 
(V, Mo) that co-ordinate with the titania support. Jibril et al. [ 27 ] studied alkali and 
molybdenum promoted Cr203/Al203 catalyst and noted that a particular composition of 
cesium-molybdenum-chromium on AI2O3 exhibited a propane conversion of 15 . 1 % and a 
propene selectivity of 64 . 5 %. For the ethane ODH reaction phosphorus was used as a 
promoter for the Cr203/Ti02 catalyst and an optimum Cr:P ratio was required to optimize 
yields. A few other studies exist that study the effect of promoters on supported 
chromium oxide catalyst [ 28 , 29 ]. 

The above studies do not, however, provide details on the effect of promoters on 
the reaction mechanism of alkane ODH. Knowledge of the ODH mechanism is essential 
to understand the functioning of the catalyst and devise methods to improve it. One of the 
methods to improve the knowledge of the alkane ODH mechanism is by analysis of the 
reaction kinetics with the help of kinetic-parameters. Although several metal oxide based 
catalytic systems have been studied for the ODH reactions, only a few studies are aimed 
at understanding the mechanism especially by analysis of the kinetic-parameters. 
Moreover, no significant studies are reported on the kinetics, reaction mechanism and 
kinetic-parameters of ODH reaction over any modified supported chromium oxide 
catalysts to explain the effect of the modifier. 

Determination of kinetic-parameters for a particular reaction mechanism is not a 
straightforward task. It requires proper design, modeling and data collection for the 
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reactor. Furthermore, determination of the kinetic-parameters requires solution of 
ordinary differential equating initial value problems such that differences between 
observed and predicted concentrations are minimized. Several optimization techniques 
have been used to obtain the kinetic-parameters. Recently, the application of genetic 
algorithms (GAs) as an optimization tool for kinetic-parameter estimation has gained 
importance since global optima are achieved and initial guesses of parameters are not 
required. 

1.1 Objective of study 

The objective of the present study is to understand the effect of modifiers on CriOa/TiOz 
for the ODH of propane. Phosphorus was chosen as a modifier since higher fields are 
achieved. A 3wt. % Cr203/Ti02 catalyst is chosen as the base catalyst to understand the 
effect of phosphorus modification on the propane ODH reaction. The 3% Cr203/Ti02 was 
chosen since only surface active chromium oxide species are present [26]. To achieve 
the objective, a series of catalysts are synthesized by loading titanium dioxide (Ti02) with 
chromium oxide and different amounts of phosphorus. The synthesized catalysts are 
initially characterized through various characterization techniques like BET, X-ray 
diffraction, electron paramagnetic resoriance and temperature programmed reduction to 
examine the possible structural changes on the catalyst surface. The unmodified and 
modified Cr203/Ti02 catalysts are then studied for the ODH of propane at different 
operating conditions to understand the effect of phosphorus modification. The kinetic- 
parameters are estimated using a genetic algorithm and related to the structural changes. 
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Finally, the effect of phosphoms modification is proposed by obtaining the structure- 
reactivity relationship. 

1.2 Thesis organization 

The thesis contains five chapters of which the current chapter is the first. The next 
chapter. Chapter two, presents a brief literature review related to the present study. In 
Chapter three, details regarding preparation of the catalyst, characterization techniques, 
experimental set-up and the methodology for kinetic-parameter estimation are provided. 
The results obtained from the characterization and reaction studies of unmodified and 
modified Cr 203 /Ti 02 catalysts on the ODH of propane are presented and discussed in 
Chapter four. Conclusions and recommendations for future work are reported in Chapter 
five. In each chapter the Tables and Figures are given at the end in numerical order. The 
references used throughout the thesis are given at the end in numerical order as they 
appear in the text. Finally, the raw data for propane ODH over the catalysts are given in 
the Appendix. 
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Chapter-2 


Literature Review 

Supported chromium oxide catalysts have been successfully used for the ODH of alkanes 
[20, 30, 31], Unfortunately, alkene selectivity obtained by this catalyst system is very 
low. Bras et al. [24] observed the requirement of redox and Lewis acidity properties in 
the catalyst for the selective oxidation of hydrocarbons to corresponding alkenes. It has 
also been proposed that an acid-base bifunctional catalyst possessing weak acid- base 
pairs would exhibit better catalyst activity and selectivity for the ODH reactions [25], 
Previous studies suggest that the introduction of an additive or modifier to the 
unsupported metal oxide catalysts is useful for better alkene selectivity [32-37], Till date, 
however, very few studies exist that discuss the effect of modifiers for supported 
chromium oxide catalysts. A brief literature review of the influence of additives on the 
structure-reactivity relationship of supported metal oxides, in general, and supported 
chromium oxide, in particular, for the ODH of alkane is given below. This literature 
review also gives a brief detail about Genetic Algorithm (GA), which is used as an 
optimization tool for finding oiit the kinetic parameters involved with the reactions. 

Ziyad et al. [38] studied the effect of phosphorus additive on Cr 203 /Ti 02 
catalyst for ODH of C 2 H 4 and observed that the presence of chromium oxide notably 
increases the activity toward the total oxidation of hydrocarbon. The addition of 
phosphorus to the system improved the global conversion and the ethylene selectivity. An 
optimum P/Cr value at which the conversion and selectivity are best was reported. 
Characterization of the catalysts was achieved by XRD, EPR and UV-visible 


7 



spectroscopy. XRD patterns revealed that the TiOa did not suffer any structural 
modification after impregnation. EPR and UV- visible spectroscopies revealed that the 
sample essentially contained Cr^^ and Cr^^ ions, which were not very efficient in mild 
oxidation. They concluded that the improvements in catalytic activity and selectivity, due 
to addition of phosphorus, can be attributed to the appearance in the samples of 
phosphates and clustered Cr^"^ ion. 

Osiysky et al. [39] studied the chromium phosphate catalyst for the oxidation of 
n-C 2 H 4 . The kinetic regularities of the n-C+Hio catalytic transformation into CrP 04 have 
been investigated. The kinetic parameters of the n-butane oxidation reaction and the 
catalytic transformation schemes were suggested. The product maximum output 
condition under a partial n-butane oxidation on the chromium phosphate catalyst was 
suggested. 

Maiti et al. [40] investigated the effect of surface phosphorus on the oxidative 
dehydrogenation of C 2 H 6 and found out that a phosphorus-enriched surface enhances the 
ODH yield of ethane to ethylene. To understand the role of phosphorus, they carried out 
the ODH reactions on a silica surface, with and without phosphorus, using the density 
functional theory code DMol3 in a periodic super-cell. The simulation study revealed that 
activation barriers for the rate-limiting steps are lowered by 10 kcal/mol in the presence 
of phosphorus. This decrease results from a transition state in which the phosphorus atom 
remains quasi-5 -valent and fourfold coordinated. 

Sidorchuk et al. [41] determined the effect of different metal phosphates on the 
partial oxidation of C 4 H 10 by comparing the catalytic activities of Al, Cr, Fe, Ti, Zr, Sn, 
and V phosphates. It was observed that except for V phosphates, the products were 
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mainly CO and CO2. The (V0)2P207 phase exhibits maximum activity and selectivity for 
maleic anhydride formation. 

Loukh et al. [ 42 ] studied the oxidative dehydrogenation of C2H6 to C2H4 on V- 
and Cr-based phosphate catalysts at 550 °C. Vanadium and chromium cations were 
introduced in zirconium hydrogen phosphates either by cationic exchange of their acidic 
proton or by impregnation with and Cr^"^ and then compared to (V0)2P207 and 
CrP 04 pure phases. They reported that and Cr^"^ on the aZr(HP 04 ) 2 . 2 H 20 phase 
were more selective toward ethylene. The (VO)2P207 and to a lesser extent CrP 04 were 
even more selective toward C2H4 at similar conversion levels. UV, ESR, XRD, IR, and 
SEM were used for the characterization of the catalysts after catalytic testing at 823 K. 
Based on the obtained data they concluded that the catalytic features were related to V 
and Cr local arrangements (small clusters or chain arrangements) and to the counter 
anion, O^' or P 04 ^’. The best catalyst for ethane oxidative dehydrogenation consisted of 
and Cr^"^ chains separated, by P04^' anions. 

Huff and Flick [ 14 ] examined the ODH of C2B6 using Cr203 and Pt-modified 
Cr203 -coated ceramic foam monoliths and Cr203/Al203 pellet catalysts. They observed 
the high C2H4 selectivity and C2H6 conversion over the Cr203 catalysts compared to Pt- 
coated monoliths. The 10 wt.% Cr203/Zr02 monolith catalyst produced the best results of 
all catalysts examined, giving C2H4 selectivities of 70 % at a C2H6 conversion of about 
80 % compared to a C2H4 selectivity of 65 % at a C2H6 conversion of 75 % under the same 
conditions over the Pt/OK-Al203 monolith catalyst. 

Ramis and Busca [ 43 ] reported the effect of several modifiers on supported 
chromium oxide catalysts. For example when an alkaline earth metal was added, the 
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lewis acidity decreased whereas non-metallic dopant (S04^', PO"^') increased the lewis 
acidity. Elements such as etc. which are non-interactive will not affect 

the structural of surface chromyl species, however, will be competitive with the basic 
sites of the support in the reaction. The effect of additives on the reactivity of each 
supported catalysts will be unique as their influence on the activity depends upon the 
reaction mechanism. 

Neri et al. [44] investigated the oxidative dehydrogenation (ODH) of iso-C 4 Hio on 
Ca-doped chromium oxide catalysts supported on AI2O3. They observed dispersed Cr^"^ 
and Cr®"^ species in addition to Cr 203 and CaCr 04 crystallites on the alumina surface. The 
relative amount of the chromium species depends on the Ca loading. The Cr^VCr^"^ ratio 
decreases on increasing the Ca loading due to the preferred formation of bulk chromate 
species. The Ca loading affects the reducibility of the Cr®"^ species and the acid sites 
strength distribution of the catalysts. 

Gr 2 ybowska et al. [11] performed ODH of iso-C 4 Hio on chromia supported on 
Si02, AI2O3, Ti02, Zr02. and MgO. They also examined the effect of potassium- 
modification. The supported catalysts have found to be active in the ODH of isobutene at 
relatively low temperatures 200-400°C and the total activity and selectivity to isobutene 
is dependent on the nature of the support and the potassium presence. The highest 
selectivity to isobutene (70% at 5% conversion) and the isobutene yields of 9% were 
obtained for CrOx/ Ti02 and K- promoted Cr 0 x/Al 203 catalysts. 

2.1 Kinetic Modeling of ODH 

In the last few years several studies were made on the structure and reactivity of ODH of 
propane over supported chromia and vanadia catalysts. However, only a few kinetics and 
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mechanistic studies on the oxidative dehydrogenation of propane over supported metal 
oxides are reported. A brief literature survey on the kinetic study for ODH and kinetic 
parameter estimation is given below. 

2.2 Kinetic Study for ODH 

Kinetic isotopic studies were carried out with supported V 2 O 5 catalysts assuming 
a first order dependency with partial pressure of propane [45]. It was observed that CO 2 is 
formed by primary decomposition of propane and secondary decomposition of propene, 
whereas CO is formed by the secondary decomposition of propene. A Mars-Van 
Krevelen redox mechanism involving two lattice oxygens in irreversible C-H bond 
activation was proposed. In these reaction mechanisms, the reaction rates were proposed 
to be independent of oxygen partial pressures. Chen et al. [46] suggested that the lattice 
oxygen in the primary dehydrogenation and combustion were different. Chen et al. [47] 
investigated ODH of propane over supported V, Mo and W oxides and concluded that the 
mechanism was the same for all these catalysts. It was suggested that the propane ODH 
and propene combustion rates depends upon the C-H bond energy difference and also on 
the adsorption enthalpies of propene and propane. 

Creaser [48] studied the partial pressure dependencies of propane oxidation. The 
corresponding rate expression, which is in accordance with Rideal type rate expression 
for various concentrations of vanadium on amorphous AIPO4, is adequate to explain the 
partial pressure dependencies of both oxygen and propane. A first order dependency of 
propane for the ODH reaction is proposed. 
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Creaser and Andersson [49] have undertaken a detailed kinetic investigation of 
ODH on V-Mg-0 catalyst. Various rate expressions are derived in order to fit the data 
based on PoAver Law type models and mechanistic models, such as, Langmuir- 
Hinshelwood and Mars-Van-Krevelen models. In these models, carbon oxides were 
considered as secondary products. 

Grabowski et al. [50] tested various kinetic models on the propane oxidative 
dehydrogenation on the VaOj/TiOa catalyst. Numerical methods (the Runge-Kutta 
method combined with the Levenberg-Marquardt method) were applied to solve the 
systems of differential equations describing the kinetic of the ODH. It was observed that 
the steady-state model, in which the surface oxygen plays an important role, provided the 
best description of the experimental results. The analysis of the reaction network 
indicated that propene is the only useful product primary produced, and that CO and CO 2 
are produced largely by the sequential oxidation of the in situ produced propene, and to a 
lesser extent by a parallel route of the direct deep oxidation of propane. 

Bottino et al. [51] studied the kinetic of propane ODH on a V/y-AllOS catalyst 
prepared by the adsorption technique. The kinetic parameters were determined through 
non-linear regression analysis. Two approaches were applied. In the first one, power law 
expressions were used to describe the network of reactions involving propane, propylene, 
CO and CO 2 . In the second approach; the experimental data was fitted with a rate 
equation obtained by assuming that molecular oxygen re-oxidizes the active site reduced 
by the surface. 

Parameter estimation using traditional methods, such as, Levenberg-Marquardt’ s 
method depend strongly on initial guesses values. On the other hand, objective functions 
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based on non-linear models are common features in most recent parameter estimation 
problems. These nonlinear models generally contain more than one minimum, 
consequently, requires a tedious algorithm, which sometime reaches to non-global 
optima. Thus, there is a need for development of efficient algorithms to replace the 
traditional methods. One such class of algorithms known as Genetic Algorithms (GA) 
fulfils this need satisfactorily and shows promising results in terms of its accuracy and 
efficiency. 

2.3 Kinetic Parameter Estimation 

Conventionally, nonlinear regression method has been used for obtaining kinetic 
parameters but its use for a case of multi-response and correlated systems is not an 
appropriate one [52]. It is proposed that a determinant criterion can be considered to be 
the best one for obtaining the Idnetic parameters for the above mentioned case. Genetic 
algorithm has recently replaced the nonlinear regression methods to optimize the 
determinant. Genetic algorithm (GA) is used to find approximate solutions to difScult- 
to-solve problems. , 

In an attempt to find initial estimates of rate constants for non-linear chemical 
kinetics. Wolf et al. [53] have exploited GA techniques without a priori assumptions of 
rate determining step in order to apply for the wide range of partial pressures of reactants 
and temperatures. Elliot et al. [54] have used inversion procedure of genetic algorithm to 
estimate rate parameters corresponding to product species measurement data from 
combustion of fuel experiments. This study suggests that its wide application to other 
chemical kinetics and optimization of other higher order kinetics is possible. Polifke et al. 
[55] employed genetic algorithm to investigate the kinetic parameters of simplified 
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reaction mechanism for methane combustion. It is also determined that using genetic 
algorithm requires minimum human effort and little insight into the details of chemical 
mechanism to generate reliable kinetic parameter. 

Several of studies have been made to implement GA for obtaining kinetic 
parameters. The performance of genetic algorithm strongly depends on the genetic 
algorithm running parameters. Furthermore, it is observed that the application of genetic 
algorithm to generate a good initial point for the non-linear local convergence method 
requires less computing time and increase the reliability of the parameters. 

2.4 Genetic algorithm (GA) 

GA algorithms mimic evolution and natural selection to solve a problem [56-58]. Due to 
this the stronger individuals are likely to be the winners in a competing environment. 
GA’s are gaining immense popularity for their potential as an optimization technique that 
can effectively search in the bounds of the complex problems. These are multidirectional 
search and optimization procedures that combine survival of the fittest among a set of 
string structures with structured yet randomized information. Optimal binary strings are 
located by processing an initially random population of strings using artificial mutation, 
crossover and selection operators, in an analogy with the process of natural selection. In 
this way better models or solutions can be bred from an originally random starting 
population. 

Initially a population of individuals (solutions) is created, each represented by a 
chromosome (a collection of genes or characteristics) to form an initial pool of possible 
solutions for the given objective function to be optimized. This is called the first 
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generation pool. All of the chromosomes, having a particular fitness value, are evaluated 
by an objective junction and effectively the pool is sorted with those having better fitness 
at the top, representing better solutions to the problem. 

The next step of the algorithm is to generate a second generation pool of 
parameter lists, which is done using the genetic operator, selection, crossover and 
mutation. 

• Selection: The first step in the construction of the next generation is to select a 
pair of chromosomes for crossover by some well-defined selection methods like roulette 
wheel selection, tournament selection etc. Tournament selection is used in the present 
study. Selection is more biased towards elements of the initial generation which have 
better fitness value. This process is repeated until the new offspring pool is full. 

• Crossover: Following selection, the crossover (or recombination) operation is 
performed upon the selected chromosomes. In general, two parent chromosomes are 
chosen to reproduce and their crossover results in two new child chromosomes, which are 
added to the second generation pool. This is repeated until there are an appropriate 
number of candidate solutions in the second generation pool. 

• Mutation: The next step is to mutate the newly developed pool, again by a 
process of selection, this time of individual chromosomes, followed by application of 
the mutation genetic operator. The main role of mutation is to provide genes not present 
in the initial population so as to prevent stagnation at local optima. 
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These processes result in a second generation pool of chromosomes that is 
different from the initial generation, which is then evaluated and the fitness values for 
each list is obtained. Generally the average degree of fitness will have increased by this 
procedure for the second generation pool. The process continues by generating third, 
fourth, fifth and future generations, until one of the generations contains solutions which 
are acceptable. 

A flowchart explaining all the major operations of GA as briefly discussed above 
is shown in Fig. 2 . 1 . 

2.5 Summary 

Thus, previous studies suggest that better selectivity and reactivity performance of 
supported chromium oxide and other metal oxides modified with additives like P, K, Ca, 
La is possible. Furthermore, phosphorus improves the activity and selectivity of the 
catalysts in the ODH of alkanes by introducing some specific redox and acid-base 
properties in the supported compound. In the present study, to understand the effect of 
phosphorus as a modifier, a P-Cr203/Ti02 with different chromium to phosphorus molar 
ratios are synthesized and then its effect on the structure and reactivity behavior of the 
catalysts is studied. By estimating the kinetic parameters also, an understanding of the 
role of phosphorus as a modifier, for the ODH of propane, is developed. 
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Figure 2.1 A flowchart of the working principle of GA 
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Chapter-3 


Experimental Details 

3.1 Sample Preparation 

The incipient wetness impregnation method was employed to synthesize the titania 
supported unmodified and modified chromium oxide catalyst samples. This method is 
knoAvn for its simplicity and viability for making low amounts of metal oxide catalysts 
[59], Modified catalysts were prepared by co-precipitation of chromium and phosphorus 
ions. Titania (TiOa, Degussa, P-25) was used as a support. Chromium nitrate 
nonahydrate (Cr(N03)3 -91120, Aldrich, 99.98% purity) and diammonium hydrogen ortho- 
phosphate ((NH 4 ) 2 HP 04 , S. D. Fine-Chem Ltd.) were used as the precursors for the 
chromium and phosphorus, respectively. Initially, the titania support was pretreated with 
incipient volumes of double distilled water. The support was then dried at room 
temperature in a desiccator followed by drying at 383 K for 8 h. Finally, the support was 
calcined at 723 K for another 8 h. 

The incipient wetness volume of solution was prepared by adding precalculated 
amount of chromium nitrate nonahydrate and diammonium hydrogen ortho-phosphate to 
double distilled water. The mixture containing chromium precursor, phosphorous 
precursor and water was stirred until a homogeneous solution was formed. For preparing 
the unmodified chromia-titania (CrTi) catalyst the phosphorus precursor was not 
required, and similarly, for the phosphorus-titania (PTi) catalyst the chromium precursor 
was not required. The above solution containing chromium and/or phosphorus ions was 
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intimately mixed with the pretreated support to form a paste. The paste was dried at room 
temperature for 12 h, followed by drying at 383 K for 12 h and finally calcined at 723 K 
for 12 h. 

In the present study five catalyst samples were prepared. The nomenclature of the 
CrTi and PTi samples were based on the weight percent of the CraOs and P 2 O 5 present in 
the sample. The nomenclature of the modified CrTi catalysts was based on the molar 
ratio of chromium to phosphorous on the TiOa support. The nomenclature and 
composition of the samples synthesized are tabulated in Table 3.1. 

3.2 Characterization 

3.2.1. Surface Area Studies 

The surface areas of the samples were obtained using the multi-point BET method. A 
bench-top COULTER SA 3100 analyzer equipped with SA-VIEW TM software using 
Na adsorption at 77K was used for this purpose. 

3.2.2. X-Ray Diffraction (XRD) Studies 

X-Ray diffraction patterns were obtained in the range of 20-80°(scanning rate: 10 ° min 
with a 180 Debye flex-2002 X -ray powder diffractometer equipped with a 
monochromator. Ni filtered Ka radiation from a Cu target (^= 1.54056 A°) was used. 

3.2.3. Temperature Programmed Reduction (TPR) Studies 

The TPR studies were carried out in a microreactor containing 0.03 5g of the sample and 
attached to a Micromeritus Pulse Chemisorb 2705 analyser. The same amount of each 
sample was used for the TPR studies. Helium was used as a carrier gas and also to degas 


19 



the samples prior to the reduction experiments. The carrier gas flow was adjusted to a 
flxed value between 30 and 50 cc/min. Degassing of the sample was achieved at 523 K 
for 0.5 h. A 5% YyAx mixture flowing at 40 seem was used for reduction and the 
temperature was ramped at a rate of 10 K/min from - 373 to 973 K. The amount of 
hydrogen consumed was detected using a TCD. 

3.2.4. Electron Paramagnetic Resonance (EPR) Studies 

The EPR spectra were obtained in aBRUKEREMX 1444 EPR spectrometer. The spectra 
were obtained under ambient conditions using a microwave frequency of ~ 9.86 GHz and 
a microwave power 0.20 mW. The scan time was 8 min and the magnetic field 
modulation frequency was 1 00 kHz. The microwave radiation travels down a waveguide 
(a type of rf pipe) to the sample, which is held in place in a microwave ‘cavity’ 
positioned between the poles of two magnets. Spectra are obtained by measuring the 
absorption of the microwave radiation while scanning the magnetic-field strength. The 
EPR spectra are displayed in derivative form to improve the signal-to-noise ratio and 
were calibrated with DPPH using a duel cell. 

3.3 Reactivity 

3.3.1 Experimental Set-up 

The ODH of propane (CsHe) is performed in a fixed bed, down-flow, tubular quartz 
reactor. The reactor was a single piece of quartz with an inlet of 10 mm inner diameter 
and an outlet of 5mm outer diameter. Schematic of the reaction set-up is shown in Fig. 
3.1. 
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As shown in Fig. 3. 1 two separate thermal Mass flow Controllers (Bronkhost Hi- 
Tech, Model F-201D and Model F-201C) were used to individually control the flow 
rate of propane and air. The two gases are mixed at a T-junction and then sent to the 
reactor. The reactor tube was mounted vertically in a tubular furnace. The catalyst bed 
containing the catalyst and quartz powder mixture was placed on quartz wool at the 
center of the 300 mm long quartz reactor. Quartz powder acts as a diluent with the 
catalyst to prevent temperature gradients and to avoid the channeling of gas within the 
catalyst bed. A chromel-alumel thermocouple was inserted from the top of the reactor to 
measure the temperature of the reactor and the catalyst bed. It was placed in such a 
marmer that the tip of the thermocouple remains just above the catalyst bed. The 
temperature of the catalyst bed was controlled by a PID temperature controller (FUJI 
Micro-Controller X Model PXZ-4). 

Product gases were analyzed using an online gas chromatograph (AEMEL 
NUCON 5765) using a Hysep-Q column for separation of hydrocarbons, CO and CO 2 . 
Reaction products were detected with the help of a Flame Ionization Detector (FID). 
Propene (CaH^), CO and CO 2 were found to be the main reaction products. A methanizer 
attached with gas chromatograph converts the carbon oxides (CO and CO 2 ) to methane in 
order to be detected by the FID. No oxygenated products were found. Several runs were 
taken at each experimental condition to ensure that steady state conditions was achieved. 

Blank experiments were also performed to determine the homogeneous 
contribution to the reaction. The fact that no reactions occurred in the blank experiments 
indicated that gas-phase-initiated reactions did not occur at the reaction conditions 
considered, hi addition, reactions in the gas phase involving desorbed radical species 
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were not considered to play a role due to the moderate reaction temperatures and 
precautions taken to minimize the void gas volume in the reactor. 

3.3.2 Reaction Studies 

Two types of reaction studies were considered: 

i. Contact time variation studies 

ii. Operating conditions variation studies. 

3.3.2. 1 Contact time variation studies 

Contact time studies for ODH of propane were performed using 0.2g of 3% Cr 203 /Ti 02 
and phosphorus modified CrTi catalyst samples at 2:1 C3H8:02 ratio with a total flow rate 
ranging from 20 to 120 cc/min. Temperature was kept constant at 653 K for all the 
catalysts. For each operating condition two runs were taken and average was reported. 

The effect of contact time, W/Fad, on the conversion and selectivity was studied 
where W is the weight of the catalyst and Fao is the total flow rate. 

3.3.2.2 Operating conditions variation studies 

The operating conditions considered were C3H8:02 molar ratio and temperature. A 
physical mixture of 0.2 g of the catalysts and required amount of quartz glass powder to 
form a bed height of 1 cm was used for this reaction study. Experimental runs were 
performed at three different temperatures of 633, 653 and 673 K with a constant total 
flow rate of 75 cc/min. The C3H8:02 flow ratios were varied from 1:1 to 3 :1. 

3.3.2.3 Analysis of reaction data 
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The products analyzed from the reaction outlet were CaHs, CO, CO 2 and unconverted 
CsHg. Each of these components gave separate areas in the chromatograph, which were 
then converted to mole fractions based on their response factor. The mole fractions were 
then used to calculate the conversion, selectivity and yield per gram catalyst as per the 
formula given below. 

i. Conversion 

Conversion for the reaction was calculated by using the following formula. 

Propane Conversion (%) =(nc/nf)xl00 (3.1) 

ii. Selectivity 

Product selectivity for the reaction was calculated by using the following formula. 

Product Selectivity (%) = (nhc / ik) x (Nhc/ N?) x 1 00 (3.2) 

iii. Propene Yield 

The amount of propene formed based on the inlet propane was represented by yield, 
which was calculated as follows: 

Propene 5 aeld (%) = (nhc / nf) X ( Nhc/ Np) X 1 00 (3.3) 

iv. Carbon Balance 

Carbon balance represents the accuracy of measurements and it provides an indicator as 
to whether carbon is deposited on the catalyst. Carbon balance calculations are done 
based on the following formula. 

Carbon balance = (Nhc/Np) (3.4) 


23 



Where, 


tic number of moles of propane consumed 
nf number of moles of propane fed 

nhc number of moles of products formed. The products formed are propene, 
C0andC02 

Nhc the number of carbon atoms present in the products 
Np the number of carbon atoms present in the propane 

3.3.3 Modeling of the Reactor and Kinetic Parameter Estimation 

Modeling of the reactor is achieved by the integral method. With appropriate reactor 
modeling the kinetic-parameters are estimated for steady state reactor operation. 

3.3.3. 1 Modeling of the reactor and problem formulation 

To model the reactor the following assumption were made. 

> The reactor operates under isothermal and steady state condition 

> Total number of moles of gases remain constant 

> There is no heat and mass transfer limitations 

> Gas phase reactions are negligible 

> Catalyst activity remains constant. (Negligible catalyst deactivation) 

All assumptions are justified since propane conversions are low and feed gases are 
diluted with nitrogen. 

The differential material balance equation for each component, i, for a particular 
reaction network can be written as 

Ky, U -Ky, ^ P.s) 
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Where, 


Vo = volumetric flow rate of the feed, mVs 
yi = mole fraction of the i component 

rii} = stoichiometric coefficient of the i'** component for the reaction 

rj = rate of j reaction and is a function of Kj or rj = :^(yi, Kj) 

i\\ 

Kj = kinetic parameters for the j reaction 
V = number of components 
W= weight of the catalyst, kg. 

The above equation with the assumptions mentioned can be simplified and written as 
follows. 


(fyj 

dW 



(3.6) 


Equation (3.6) is a set of ordinary differential equations initial value problems (ODE- 
rVPs). The set of ordinary differential equations represented by Eqn. (3.6) can be 
simultaneously solved to obtain the output mole fraction of each component based on the 
input initial value of mole fraction of each component, yio, and temperature, provided the 
weight, stoichiometric number, n,/, volumetric flow rate, Vo, the function dependence q = 
fj (yi, Kj) and kinetic parameters are known. The input mole fraction is changed during 
integration over the entire mass using Runga-Kutta fourth order technique to give the 
output mole fraction. 

The rate of reaction, rj, is a function of kinetic-parameter, Kj. For a particular set 
of input kinetic-parameters, the output mole fraction can be determined by integrating the 
above equation. Thus, comparison of the actual output mole fraction, yi,expt, with those 
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predicted by assuming certain values of Kj, yi,predicted, can be used as a basis to calculate 
Kj. Since the kinetic-parameters have Arrhenius dependence on temperature non-linear 
regression analysis is required to obtain rate constants. 

The response data, the output mole fraction of CsHg, C 3 H 6 , CO and CO 2 , are 
assumed to be well described by a nonlinear model given by 

Yh=gh(6=y)+2:h, h=1.2.....,v (3.7) 

where, 

Yh = vector of random variables representing the response variable 
Zh = error associated while calculating the response h 
0 = parameter vector 
y = concentration vector 
gh(®>y)“ predicted concentration 

The predicted output mole fraction of the i* component in the u*** experiment, 
yiu,pred obtained this way can then be compared with the actual output mole fraction 
yiu, actual- An objective function involving the predicted and actual mole fractions of the 
components is minimized to obtain the best value of the kinetic parameters. 

For multiresponse systems and when responses are correlated, the ideal criterion 
is the minimization of the following determinant [60-63]. 
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^ (yiu ” yiu, pred) 
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The problem then can be formulated as minimization of the above determinant. 
Minimization of the deteraiinant is achieved by applying Genetic Algorithm. 

3.3.3.2 GA - Optimization Technique 

The optimization technique applied in the present study is a Genetic Algorithm (GA). It 
is considered to be robust and efficient method for optimization. The performance of 
genetic algorithm strongly depends on the genetic algorithm running parameters as 
discussed in Chapter-2. Here a brief detail of source code and the GA parameters used in 
the present study to estimate the kinetic-parameters is given. 

3.3.3.2.1 Source Code 

A NSGA code (Non-dominated Sorting Genetic Algorithm) has been used in the present 
study to determine the kinetic parameters. The NSGA code is taken from the KANGAL 
lab of nT Kanpur, and has been developed by Deb and co-workers. 

The GA implementation is restricted to real coded variables only. All constraints 
are greater-than-equal-to type (g >= 0) and normalized. Evaluation is done based on the 
determinant criterion. The above determinant is coded and the best values of the kinetic 
parameters are estimated by minimization of it. 
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3.3.3.2.2. GA Parameters 


Initially a population of solution is randomly generated. The size of population was 
chosen in between 120 and 200 for power law type models. The variable boundaries are 
fixed for each parameter and fixation of boundaries is done according to previous 
knowledge of the parameters. Selection of parents from the pool of solutions is done by 
tournament selection. The size of tournament in the present study was 2. Crossover 
probability exchanges information among parent solutions, whereas polynomial mutation 
operator is utilized to introduce the extra diversity to the solutions. The values of 
crossover and mutation probability were 0.8-0. 9 and 0. 1-0.2, respectively. Simulated 
binary crossover operator (SBX) is used; In this real coded GA that chases the function 
optima wherever it goes in the search space. The exponents of SBX and mutation were 2 
and 200 respectively. The termination criteria used in this GA code was the number of 
generation. Number of generations was varied from 1,000 to a maximum of 30,000. 

3.4 Reaction Models 

The kinetic study of propane ODH involves various proposed reaction models. These 
proposed models can be broadly classified as Power law type and mechanistic type. 
According to these models propene is a primary product of propane. However, the 
formation of carbon oxides as primary, secondary or tertiary products is still debatable. 
Here discussion is confined to Power Law based models and details about mechanistic 
models, such as Langmuir-Hinshelwood model, Rideal type model and Mars-van- 
Krevelen model are not considered. 
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3.5 Reaction Scheme 


The generalized reaction network for the ODH of propane is shown in Fig.3.2. 
Here the reaction of the carbon containing molecules with oxygen is presented in the 
form of six reaction steps, r; to re. Propene is formed by the ODH of propane (ri). 
Propene can degrade to form CO and CO 2 by reactions and rs, respectively. The carbon 
oxides, CO and CO 2 , can also be formed directly from propane by reactions r 4 and rs, 
respectively. Finally, CO 2 can also be formed by the oxidation of CO by reaction te. 

The stoichiometric equations corresponding to each reaction step in the reaction 
network is given by equations (3 .9) to (3.14). 


C,H,+0.5O,-^C,H,+H, 

(3.9) 

C,H,+30, -^3CO + 2H^O 

(3.10) 

+ 4 . 5 O 2 -> 3 CO 2 +3H^O 

(3,11) 

+50^ ^3CO^ +4H^O 

(3.12) 

C,H,+3.5 0^^3C0 + 4H^0 

(3.13) 

(7(9 + 0.502 — y CO 2 

(3.14) 


3.6 Power Law Models 

In the power law models it is assumed that the reaction rate is proportional to the 
partial pressure of reactant raised to a particular exponent. Based on the reaction steps 
given by Eqns. 3 .9 to 3 . 14 the reaction rates are given by Eqns. 3 . 1 5 to 3 .20 . 


29 




(3.15) 


(3.16) 


(3.17) 

=K(Pc,«.r(Po,y' 

(3.18) 


(3.19) 

=kAPcoy(Po,)’“ 

(3.20) 


Where, 

ki is the it** rate constant for reaction i and is a function of the pre-exponential 
factor, kio, and activation energy, £, 

ai and hi are the partial pressure exponent of the reactants for reaction i. 

Thus, each reaction, r,, involves four parameters; ko, Ei, ai and bi. 

Based on previous studies only two power law models were analyzed to determine the 
empirical reaction rates. These two power law models PL-1 and PL-2 are described 
below. 

3.6.1. PL-1 Model 

The PL-1 model assumes that CO and CO 2 are secondary products formed from propene. 
Thus, ri, 12 and rs are the only reaction steps considered. There are 12 parameters in this 
consecutive reaction model; four for each reaction step. This type of reaction scheme has 
been considered previously [49]. A schematic of the PL-1 model is given in Fig. 3.3. 

3.6.2. PL-2 Model 
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The PL-2 model takes into consideration ri, r 4 and rs reactions. Thus, CsHs reacts with 
oxygen and produces CsHk, CO and CO 2 from ri, r 4 and rs, respectively. Consequently, 
C 3 H 6 , CO and CO 2 are primary products. There are 12 parameters also in the PL-2 model. 
A schematic of the PL-2 model is given in Fig. 3.4. The parameters involved in the PL-1 
and PL-2 models are given in Table 3.2. 

/ 

3.6.3. Reparameterisation of reaction rates 

Parameter estimates obtained by fitting models to data are often highly correlated with 
each other. Under such circumstances the parameters are of little use for predicting the 
nature of the system. In chemical kinetics high correlations are frequently encountered 
between estimates of the kinetic-parameters in the Arrhenius expression for a rate 
constant making the analysis of a reaction network very difficult [64]. For both models 
considered here the reparameterized rate constants are expressed in Arrhenius form as 

=^,oexp[-^(^-^)] (3.21) 

m 

Where, 

kio = pre-exponential factor 
Ei =activation energy for the reaction i, kJ/mol 
R = universal gas constant, kJ/kmol K 
T= actual reaction temperature, K 
Tm = mean temperature, K 
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The inverse temperature is centered about the inverse mean temperature (1/Tm). 
Furthermore, the partial pressure of a component, Pj, is also centered about the mean 

P. 

partial pressure of the component, i.e.,(-^)‘'' . 


Where, 

P" = mean partial pressure of the component j, atm 

This centering procedure eases the parameter search by minimizing the 
statistical correlation between activation energy and pre-exponential factor [64, 65], 
Thus, the reparameterized rate equation is given by 


a;. =A:,,exp[-^(^--;^)] 


Ml 

^1 

( p 

J 


p ^ 

\^k J 


(3.22) 
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Table 3.1 Composition and Nomenclature of TiOj supported Cr-P-O catalysts 


Nomenclature 

Cr : P molar 

ratio 

% CrzOa 

loading 

%P205 

loading 

T 1 O 2 

- 

- 

0 

3CrTi 

- 

3 

0 

2CrlPTi 

2;1 

2.96 

1.38 

l.SCrlPTi 

1.5;1 

2.95 

1.83 

ICrlPTi 

1:1 

2.92 

2.72 

1.4PTi 

- 

- 

1.4 
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Table 3.2 Different Power Law models and corresponding number of kinetic 
parameters 


Kinetic 

Contributing 

Number of Parameters 

Model 

reactions (Fig. 3.2) 

kio 

Ei 

ai 

bi 

Total 

PL-1 

ri,r2,r3. 

3 

3 

3 

3 


PL-2 

n, r4, rs 

3 ! 

3 

3 
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Figure 3.2 Generalized reaction network: Carbon containing 
molecules only shown 
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Fig 3,3 PL-1 reaction model 



Fig 3.4 PL-2 reaction model 
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Chapter-4 


Results and Discussion 

Several unmodified and modified Cr203/Ti02 catalysts were prepared. The unmodified 
catalysts contained 3 wt % of Cr203 on titania. The modified catalysts possessed ~ 3 wt 
% Cr203 and three different chromium to phosphorus molar ratios: 2:1, 1:5:1 and 1:1. In 
addition to these four samples, 3CrTi, 2CrlPTi, l.SCrlPTi and ICrlPTi, another sample 
containing 1 .4 vv1;% of P2O5 was also prepared. The BET surface areas of the prepared 
samples were initially obtained. The samples were further studied by XRD, EPR and 
TPR characterization techniques. The samples were then studied for the ODH of propane 
with the primary objective to understand the effect of phosphorus on the kinetic 
parameters. The results of these studies are summarized below. 

4.1 Characterization Studies 

4.1.1 Surface Area 

The surface areas of the Ti02 support, modified and unmodified chromia over titania 
samples were obtained and are tabulated in the second column of Table-4.1. Also 
included in Table 4.1 are results from TPR studies that are discussed later. It is observed 
that the surface areas were relatively constant (41 to 44 m^/g) for modified and 
unmodified Cr203/Ti02 samples and are slightly lower than the pre-treated pure support, 
46 m^/g. The decrease in surface area of samples may be due to the initial plugging of 
pores as suggested previously [66, 67]. Interestingly, when chromium oxide is deposited 
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on titania it gains stability against further decrease in surface area when the Cr203/Ti02 
sample is modified with phosphorus. Similar results were reported for vanadium oxide 
modified CxiOzIZxOi catalysts [68]. 

4.1.2 X-Ray Diffraction (XRD) 

The XRD patterns of the titania support and supported samples are shown in 
Fig.l. For all the samples peaks are observed at 20 values of 25.51, 27.61, 36.31, 38.16, 
41.71, 48.26, 54.46, 55.26, 63.16, 69.36, 70.76, 75.71 and 76.31°. These peaks are due to 
the Ti02 support. Consequently, other phases of chromium, phosphorus and/or titania are 
not detected by XRD. Small crystals of less than 40 nm size may still be present since 
they can not be detected by XRD. Similar conclusions are reported in the literature [69, 
70]. 

4.1.3 Temperature Programmed Reduction (TPR) 

The TPR profiles for unmodified and phosphorus modified chromia/titania (CrTi) along 
with phosphorus/titania (PTi) samples are shown in Fig. 4.2 where the TCD signal is 
plotted versus the temperature. It is observed from Fig. 4.2 that for all the phosphorus 
modified and unmodified CrTi samples two Tmax temperatures were observed. The first 
reduction peak is observed at ~615 K and the second peak is observed to shift as the 
phosphorus content is increased. For 3CrTi the second reduction peak is at 736 K that 
shifts to 774 K for 2CrlPTi, 809 K for 1.5CrlPTi and 860 K for ICrlPTi. A high 
temperature reduction peak at 890 K is detected for 1.4PTi. The reduction peak at -615 
K, present in all chromia containing samples, corresponds to the reduction of the surface 
chromium oxide species [22]. Two reduction peaks were observed for Cr203/Al203 and 
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Cr 203 /Si 02 samples also [71]. The Tmax values for the samples are given in the third 
column of Table 4.1. 

For quantitative analysis, the peak areas were integrated and the amount of 
hydrogen consumed was determined by injecting known amounts of hydrogen. From the 
calculated peak area of the first reduction peak and the calibration amounts of hydrogen, 
the hydrogen consumption is calculated. Based on the hydrogen consumed and the 
amount of chromia present the H/Cr molar ratio is calculated. The H/Cr ratios are 
tabulated in the fourth column of Table 4.1 for the unmodified and modified CrTi 
catalysts. From the data presented in Table 4.1, it appears that the first Tmax temperature 
appears to be relatively constant. Furthermore, the H/Cr ratio appears to decrease with the 
addition of phosphorus to the CrTi catalyst suggesting that the amount of the surface 
chromium oxide species decreases. 

4.1.4 Electron Paramagnetic Resonance (EPR) 

The EPR spectra are obtained under ambient conditions for modified and unmodified 
chromia/titania samples and are shown in Figure 4.3. EPR spectroscopy is known to be 
sensitive to Cr^"*" (d’) and Cr^"^ (d^) cations and not to Cr^"^ (d°, diamagnetic) ions [72-75]. 
The catalyst free of phosphorus, i. e. the 3CrTi sample, exhibits typical strong axially 
symmetric peak with the parameters g=l. 98 and AH= 50 G. This peak corresponds to the 
y-signal of isolated axially symmetric Cr^"^ species [76]. The Cr^"^ species is also found in 
the EPR spectra of the phosphorus modified CrTi samples. On adding phosphorus to the 
CrTi sample, the y-signal intensity decreases and a broad p-signal with g = 1.98 and A H 
= 820 G appears. This broad signal was also observed by Loukah et al. in characterizing 
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the CrP 04 phase and was attributed to the presence of clustered Cr^'^-O-Cr^'^ species in 
amorphous CrP 04 [42], The intensity of the P-signal increases with phosphorus addition. 
Consequently, with phosphorus addition the amount of Cr^"^ species decreases and the 
amount of Cr^'^-O-Cr^'^ species associated with amorphous CrP 04 phase increases. No 
peak was found in the spectrum of 1.4PTi and is not shown. 

Thus, characterization studies reveal that while preparing the catalyst the support 
is not significantly affected. Surface reducible chromia species are present in the 
modified and unmodified CrTi samples in addition to a Cr^"^ and Cr^"^ species associated 
with amorphous CrP 04 . The amount of surface reducible chromia species decreases and 
the Cr^"" species increases with an increase in phosphorus content. 

4.2 ODH Reaction Studies 

The propane ODH reaction with the unmodified and modified CrTi catalysts was 
performed. The unmodified and phosphorous-modified CrTi samples were observed to be 
active for the propane ODH reaction. Propene and carbon oxides (CO and CO 2 ) were the 
only reaction products detected. Titania and titania supported phosphorous were, 
however, inactive under the experimental conditions considered in the present study. 

4.2.1. Contact Time Study 

Initially, the propane ODH reaction was performed at a constant temperature of 653 K 
and a CsHg'.Oi molar flow ratio of 2:1. The flow rate was varied from 20 to 120 cc/min. 
The peak areas of CgHg, CgHg, CO 2 and CO obtained from the chromatographs are 
tabulated in Appendix- 1. In each table there are six columns: the first column 
corresponds to the total flow rate in cc/min, the second column represents the contact 
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time (W/Fao) in Kg-m'^s and the next four columns correspond to the peak area 
corresponding of CsHg, C 3 H 6 , CO 2 and CO, respectively. 

Based on the areas the mole fractions of the components were determined from 
which the conversion and selectivity were calculated based on Eq. (3.1) and (3.2). The 
conversion and selectivity are used to generate Figs. 4.4. and 4.5. The propane conversion 
is plotted versus the contact time in Fig. 4.4. It is observed that for all of the catalysts the 
propane conversion increases with an increase in contact time. Furthermore, at a constant 
contact time the conversion decreases as the phosphorus amount in the catalyst increases. 
Thus, the activity follows the order 3CrTi>2CrlPTi>1.5CrlPTi>lCrlPTi. 

In Fig. 4.5 the propene selectivity is plotted against propane conversion. From 
Fig. 4.5 it is evident that the propane conversion and propene selectivity are inversely 
related. Furthermore, at iso-conversions the C 3 H 6 selectivity increases with the amount of 
phosphorus in the catalyst except for the ICrlPTi sample where the selectivity decreases. 
Consequently, at iso-conversion the propene selectivity follows the order 
1.5CrlPTi> ICrlPTi >2CrlPTi>3CrTi. Thus, it appears that there exists an optimum CnP 
ratio where the selectivity, and, therefore the yield is maximum. 

4.2.2. Reaction Data for Kinetic Investigation 

The propane ODH was carried out over the active catalysts by varying the C3H8:02 molar 
ratio and temperature. The C3H8;02 ratio was varied as 3:1, 2:1 and 1:1 eind the 
temperature was varied as 633, 653 and 673 K. During these reaction studies the total 
feed flow rate was maintained at 75 cc/min. All the runs were carried out in the 
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temperature range of 633 to 673 K. Each run was repeated twice to check the consistency 
of the result. 

The raw peak areas of all the four analyzed components (CO, CO 2 , C 3 H 6 and 
unreacted CsHg) obtained during reaction over 3CrTi, 2CrlPTi, l.SCrlPTi and ICrlPTi 
catalysts are presented in Appendix-2 as Tables A-2. 1 to A-2.4, a separate table for each 
catalyst. In each table the first column corresponds to the CsHgiOj molar ratio. The 
second column of table provides the inlet area of propane, whereas the third column 
provides the temperature. The forth column of table represents the peak areas of the four 
analyzed components as obtained from the GC. The last column contains the carbon 
balance. 

The peak areas shown in Appendix-2 are used to calculate the product (CO, CO 2 
and CaHe) yields based on Eq. 3.3. The calculated yield values at all reaction conditions 
are presented in Table 4.2 to 4.5 for the four catalysts. The first column of the tables 
corresponds to the reaction temperature and the second column corresponds to the 
C 3 H 8:02 molar ratio. The product yields obtained are given in the third column of the 
table and the last column gives the carbon balance (C-balance) in terms of the ratio of 
carbon atoms out to carbon atoms in, Cout/Cin. As mentioned previously the C-balance 
provides a measure of the accuracy of GC analysis. Examination of the C-balance values 
reveals that the GC- analysis is reliable. Analysis of the yield values of Table 4.2 to 4.5 
reveals that for all the catalysts the C 3 H 6 , CO 2 and CO yields increase with an increase in 
temperature. However, at a particular temperature the yields are decreasing with an 
increase in C 3 H 8;02 molar ratio. Comparison of the yield values of one catalyst with 
another is not possible since the weight of the catalysts, and, consequently, the contact 
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times and conversions are different. Comparison of the yields at iso-conversion or the 
same contact time is more appropriate and is done based on the kinetic parameters 
estimated. 

4.3. Kinetic Parameter Estimation 

Based on the inlet CsHg and O 2 concentrations and outlet areas detected by GC 
(Appendix-2), the input and output mole percents of different components are calculated 
and shown in a tabular form in Appendix-3. There are four tables in Appendix-3 from A- 
3.1 to A-3.4 corresponding to the four catalysts. The first column in the tables 
corresponds to the reaction temperature, while the second and third columns give the 
input (C 3 H 8 and O 2 ) and output (CsHg, C 3 H 6 . CO 2 and CO) mole percents. The data in 
. each table is organized such that the C 3 H 8:02 molar ratio of 1 : 1 is the top set of the data, 
and is separated by a bold line from the following C 3 H 8:02 ratio of 2: 1 data, which is 
followed by the C3H8:02 ratio of 3:1 data. In each data set, corresponding to particular 
C3H8:02 molar ratio, the data is presented with increasing reaction temperature from 633 
to 673 K. 

Based on the input and output mole percents of CgHg, C 3 H 6 , CO 2 and CO, the 
kinetic parameters for the four catalysts based on PL-1 and PL-2 models, discussed in 
section 3.6, are determined by minimizing the objective function. As mentioned before in 
the section 3.3.2 a determinant was used as the objective function. 

4.3.1. Model Selection 

After estimating the parameters comparison of the PL-1 and PL-2 is done. Since 
the number of parameters used in the PL-1 and PL-2 models are the same comparison of 
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the determinant values used as the optimization criteria would provide the better 
suitability of the model to explain the reaction data. The determinant value (D. V.) for all 
the four catalysts obtained from PL-1 and PL-2 is presented in Table 4.6. Invariably the 
D. V. is smaller for the PL-1 than the PL-2 model for all the catalysts. Consequently, it 
appears that the PL-1 model explains the experimental data better than the PL-2 model. 
The PL-1 model suggests that a consecutive reaction takes place where carbon oxides are 
secondary products and propene is the only primary product. Similar conclusions have 
been previously proposed [49]. 

4.3.2. Kinetic parameters for the PL-1 model 

The kinetic-parameters of the PL-1 model, consisiting of pre-exponential factors, kto, 
activation energy, Ei, and partial pressure exponents, ai & bi, are determined for the 
3CrTi, 2CrlPTi, l.SCrlPTi and ICrlPTi catalysts and are presented in Table 4.7. For 
determining the kinetic parameters the values of P„^c,h^)~ 0-2958 atm, 0-148 atm 

and Tm — 653.16 K are considered. The units of the parameters are also presented in the 
tables. The last row of Table 4.7 corresponds to the P„^c,h^)= which is different for each 

catalyst. Despite PL-2 being less suitable for the present study, kinetic parameters are 
presented in Appendix-4. 

Using the kinetic parameters of PL-1 model and the experimental conditions used, 
the outlet concentration of C 3 H 8 , C 3 H 6 , CO 2 and CO can be predicted. To show the 
proper estimation of kinetic-parameters, the predicted outlet concentrations of C 3 H 6 , CO 2 
and CO are compared with the actual outlet concentrations in Fig. 4.6. From the figure it 
is clear that the kinetic parameters estimated for the PL-1 model is effective in 
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representing the propane ODH reaction under the experimental conditions used in the 
present study. 

The effect of phosphorus addition on Cr 203 /Ti 02 is analyzed by comparing the 
kinetic-parameters of the four catalysts. By examining the A:,,, values for the unmodified 

and modified catalysts in Table 4.7 reveals that the magnitude of the values decreases 
from 3.74 for 3CrTi to 1.68 for ICrlPTI as phosphorus is added. The k^o also decreases 
monotonically with phosphorus addition. The kjo value decreases for the 3CrTi, 2CrlPTi 
and l.SCrlPTi catalysts as phosphorus is added and increases for the ICrlPTi catalyst. 
For each of the catalysts k,o > kj,, > k^Q except for 3CrTi where k,Q<k 3 o>k 2 o. Creaser 
and Anderson have also observed similar trend in k,o , kjo and kjo for the propane ODH 
reaction over a VMgO catalyst [49]. 

To understand the variation of k/o values with phosphorus addition, the kjo values 
are normalized with the kw values of the 3CrTi catalyst and are plotted in Fig. 4.7 for the 
four catalysts. The normalized kio decreases with the addition of phosphorus except for 
the normalized kso value, which initially decreases and then increases for the ICrlPTi 
catalyst. However, the relative decrease of the kw values is not the same for all the 
catalysts. The normalized kw value decreases and appears to approach a relatively 
constant value. The normalized k 2 o value decreases more rapidly and almost linearly 
with phosphorus addition. The normalized kso value decreases the most rapidly and 
finally increases for the ICrlPTi sample. Thus, the pre-exponential factors do not 
decrease at the same rate with phosphorus addition. Consequently, the pre-exponential 
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factors for the propane dehydrogenation reaction and propane combustion reaction are 
affected differently with phosphorus addition. 

The pre-exponential factors are specific to the catalytic active site and depend on 
the activity per site or the number of active sites. Relating this information with 
characterization studies suggest that the number of active sites is affected since the H/Cr 
ratio associated with the surface chromium oxide species and the EPR spectra also 
changes. The characterization studies suggest that the surface chromia species are less 
and a new Cr^"^ species associated with the CrP 04 related phase is present. The decrease 
in the kio values with phosphorus addition appears to be related to the decrease in the 
surface chromium oxide species. The increase in the kio value for the ICrlPTi catalyst 
appears to be related to the formation of a new chromium oxide phase. 

The ranges of activation energy given in Table-4.7 for propane ODH, Ei, are from 
123 to 165 kJ/mol; for propene combustion to CO, E 2 , are from 24 to 55 kJ/mol; for 
propene combustion to CO 2 , E 3 , are from 22 to 61 kJ/mol. The relative trends of the 
activation energies values are consistent with previous studies that suggested E 2 and E 3 
should be less than Ej [47, 49, 51]. The difference in activation energies for ODH and 
combustion are important for obtaining high propene yields. A plot of Ei-(E 2 +E 3 } for the 
four catalysts is shown in Fig. 4.8. It is observed from the figure that the difference in 
activation energy for propane ODH and propene degradation is relatively constant for 
3CrTi, 2CrlPTi and 1.5CrlPTi, but significantly different for the ICrlPTi sample, 
consistent with the formation of new phase. 

The partial pressure exponents, ai and bi, for all the catalysts are also given in the 
Table 4.7. A pseudo-first order dependency of propane for r/ is observed for 2CrlPTi and 
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l.SCrlPTi, whereas r/ appears to be independent of propane partial pressures for 3CrTi 
and ICrlPTi samples. The oxygen dependency for the primary reaction ri is first order 
for most catalyst. The rate of CO formation is independent of oxygen partial pressure 
except for the case of 2CrlPTi where it follows a first order dependency. No specific 
trend is observed for the other exponent values. In general, it appears that there is no 
common reactant partial-pressure dependency for all the catalysts. Each catalyst has its 
specific partial pressure dependency. 

Thus, the kinetic parameters of the reactions involved in both PL-1 and PL-2 
models were successfully determined using GA for the unmodified and phosphorus- 
modified chromia over titania catalysts. The most suited model suggests that propene is a 
primary product, which then degrades to carbon oxides. Phosphorus addition primarily 
affects the pre-exponential factors, whereas the activation energies and partial pressure 
exponents are not affected or affected randomly. 

4.4. Predicted Yield 

Based on the estimated kinetic-parameters the effect of conversion on propene yield and 
selectivity can be determined. Using the kinetic parameters and for reaction conditions of 
C 3 H 8:02 molar ratio of 2; 1 the calculated propene yield versus contact time at reaction 
temperatures of 633, 653, 663 and 673 K are plotted in Fig. 4.9 for the 3CrTi catalyst. 
From Fig. 4.9 it is observed that a specific contact time exists for maximum propene 
yield. Furthermore, as the temperature increases the propene yield also increases. With an 
increase in temperature, the ratio of the rate constants for propene degradation to propane 
ODH, {k 2 + hyki, also decreases. Similar behavior is also observed for first order 
consecutive reactions where the optimum yield of the, intermediate product increases as 
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the ratio of the rate constants of the first and second reaction decreases [77]. In the 
present study the change of the ratio between the rate constants can be brought about a 
change in temperature and change in catalyst composition. 

In Fig. 4.10 the calculated propene yield is plotted versus conversion at 653K and 
again an optimum is observed except for the 2CrlPTi catalyst. For the 2CrlPTi catalyst 
the propene yield monotonically increases with conversion. Furthermore, from Fig. 4.10 
the optimum propene yield increases as the {k 2 +k 3 )lki ratio decreases except for the 
2CrlPTi where an optimum propene yield is not observed. In comparison to Fig. 4.9, in 
Fig. 4.10 the decrease in ratio is achieved by changing the catalyst 

composition. Addition of phosphorus to the CrTi catalyst decreases the ratio to 

a minimum value for l.SCrlPTi. The (Ai+A^yA/ increases again for the ICrlPTi catalyst. 

The selectivity is plotted versus conversion in Fig. 4.1 1 for all the four catalysts. 
The predicted selectivity versus conversion curves reveals that an inverse relationship 
exists between the conversion and selectivity. Furthermore, at iso-conversion the propene 
selectivity increases with phosphorus addition except for ICrlPTi sample where the 
selectivity decreases. Similar trends are observed in contact time studies in Fig. 4.5 also 
where the propene selectivity was maximum for l.SCrlPTi catalyst. 

Thus, the propene yield can be controlled to some extent by the reaction 
temperature and catalyst composition. With an increase in temperature the ratio of the 
rate constant associated with degradation to the rate constant associated with ODH 
decreases. The decrease in rate constant ratios results in a higher propene yield. An 
increase in propene yield can also be achieved by phosphorus addition since the rate 
constants ratio of degradation to ODH also decreases till an optimum Cr:P ratio of 1.5:1. 
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The decrease in propene yield for Cr:P ratio less than 1.5:1 appears to be related to the 
increase in the pre-exponential value associated with propene degradation to CO 2 . 
Interestingly, for CzHg ODH an optimum ratio of P;Cr of 1.6 was found by Ziyad et al. 
[38]. Thus, the optimum Cr;P ratio appears to depend on the ODH of the specific alkane. 
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TabIe-4.1: Surface area (S. A.), Tmax and H/Cr ratio of the catalysts 


Catalysts 

S A., 
(mVg) 

Tmax 

(K) 

H/Cr 

(molar ratio) 

Ti02 

46 

- 

- 

3CrTi 

42 

612, 736 

1.21 

2CrlPTi 

44 

613,774 

0.96 

l.SCrlPTi 

44 

613,809 

0.81 

iCrlPTi 

43 

619, 860 

0.44 

1.4PTi 

41 

890 

- 
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Table 4.2: Product yield obtained with 3CrTi catalyst for propane ODH 
Weight of the catalyst = 0.20 g; Total flow rate =75cc/min. 


Reaction 

Temperature 

(K) 

C 3 H 8 /O 2 

ratio 

Yield (%) 

C-balance 

(Cout/Cin) 

C 3 H 6 

CO 2 

CO 







633 

1:1 

0.49 

0.46 

0.15 

0.990 

2:1 

0.30 

0.44 

0.10 

0.996 

3:1 

0.22 

0.30 

0.08 

0.996 







653 

1:1 

1.40 

1.11 

039 

0.956 

2:1 

0.72 

0.8 

0.23 

0.998 

3:1 

0.48 

0.42 

0.13 

0.993 







673 

1:1 

3.97 

2.31 

1.29 

0.917 

2:1 

1.29 

0.91 

0.38 

0.989 

3:1 

1.21 

0.71 

0.37 

0.974 
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Table 4.3: Product yield obtained with 2CrlPTi catalyst for propane ODH 
Weight of the catalyst = 0.20 g; Total flow rate = 75cc/inin. 


Reaction 

Temperature 

(K) 

C 3 H 8 /O 2 

ratio 

Yield (%) 

C-balance 

(Cout/Cin) 

C 3 H 6 

CO 2 

CO 







633 

1:1 

0.26 

0.15 

0.07 

0.993 

2:1 

0.23 

0.10 

0.06 

1.002 

3:1 

0.21 

0.07 

0.05 

0.999 







653 

1:1 

0.62 

0.28 

0.13 

1.014 

2:1 

0.45 

0.27 

0.15 

0.999 

3:1 

0.43 

0.17 

0.11 

0.994 







673 

1:1 

1.33 

0.40 

0.22 

1.013 

2:1 

1.09 

0.49 

0.31 

0.989 

3:1 

0.76 

0.34 

0.24 

0.996 
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Table 4.4: Product yield obtained with l.SCrlPTi catalyst for propane ODH 
Weight of the catalyst = 0.20 g; Total flow rate =75cc/min. 


Reaction 

Temperature 

(K) 

CsHg/Oz 

ratio 

Yield (%) 

C-balance 

(Cout/Cln) 

CsHfi 

CO 2 

CO 







633 

1:1 

0.26 

0.10 

0.05 

0.998 

2:1 

0.23 

0.07 

0.05 

0.984 

3:1 

0.19 

0.07 

0.05 

0.998 







653 





0.987 




0.08 

0.973 

3:1 

0.36 

0.14 

0.09 

0.979 







673 

1:1 




0.828 

2:1 



w^m 

0.973 

mm 



mol 

0.982 




^■1 
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Table 4.5: Product yield obtained with ICrlPTi catalyst for propane ODH 
Weight of the catalyst = 0.20 g; Total flow rate =75cc/min. 


Reaction 

Temperature 

(K) 

C 3 WO 2 

ratio 

Yield (%) 

C-balance 

(Cout/Cin) 

CaHfi 

CO 2 

CO 







633 

1:1 

0.18 

0.08 

0.02 

0.992 

2:1 

0.12 

0.02 

0.02 

0.953 

3:1 

0.07 

0.01 

0.01 

0.991 







653 

1:1 

0.36 

0.15 

0.03 

0.992 



0.03 


0.944 


^5 










673 

1:1 

0.49 

0.20 

0.04 

0.996 

2:1 

0.41 

0.07 

0.06 

0.922 

3:1 

0.26 

0.04 

0.03 

0.922 















Table 4.6: Determinant values obtained during kinetic parameter 
estimation from PL-1 and PL-2 Model 


Catalyst 

Determinant Value 

PL-1 

PL-2 

3CrTi 

5.69e-27 

3.3e-26 

2CrlPTi 

3.05e-31 

9.8e-29 

l.SCrlPTi 

4.06e-29 

6.6e-28 

ICrlPTi 

5.92e-28 

5.9e-27 
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Table 4.7: Kinetic parameters for PL-1 model 


Calculated at = 0.2958 atm, 0.148 atm, T„ = 653.16 K 




Catalyst 

Units 

3CrTi 

ICrlPTi 

l.SCrlPTi 

ICrlPTi 


-o 

.f 

1 

s § 

3.74 

2.75 


1.64 

k 

^ 20 

0.77 

0.51 

0.30 

0.10 

k 

^ 30 

4.12 

1.35 

n 

1.37 

El 

O 

s 

156 

145 

165 

123 

E2 

45 

24 

55 

52 

Es 

26 

22 

40 

61 

ai 



Dimensionless 

0.18 

1.0 

1.0 

0.01 

bi 

1.0 

0.83 

1.0 

1.0 

ai 

1.0 

0.01 

0.58 

0.35 

b2 

0.01 

1.0 

0.01 

0.01 

Oi 

0,39 

0.79 

0.09 

0.58 

bs 

0.57 

1.0 

0.01 

0.55 


atm 

0.0015 

9.9E-04 

5.2E-04 

2.6E-04 
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Figure 4.1: X-Ray diffractograms of Ti02, 3CrTi, phosphorus 
modified CrTi and 1.4PTi samples 
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Figure 4.2: TPR profile for unmodified 3CrTi, phosphorus 
modified CrTi and 1.4PTi 
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Figure 4.3: EPR spectra of catalyst samples under ambient conditions 
Legends: 1- 3CrTi 2- 2CrlPTi 

3- l.SCrlPTi 4- ICrlPTi 
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Figure 4.4: Propane conversion versus contact time for unmodified 
and phosphorus modified CrTi 
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Predicted Concenteration 



Figure 4.6: Actual concentration versus predicted concentration using 
PL-1 model for unmodified and phosphorus modified CrTi 
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Catalyst, i 

Figure 4.7: Normalized k\o for all the catalysts 

Legends, i (on x-axis): l-3CrTi 2-2CrlPTi 


3-1.5CrlPTi 4-lCrlPTi 





120 



12 3 4 

Catalyst 

Fig. 4.8: Difference in activation energies of propane ODH and 
propene degradation 

Legends, (on x-axis): l-3CrTi 2-2CrlPTi 

3-1.5CrlPTi 4-lCrlPTi 
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Figure 4.9: Predicted propene yield versus contact time 


for 3CrTi catalyst; wt. of catalyst=0.2g, CaHg : O 





Propene Yield(%) 



Figure-4.10: Predicted propene yield versus propane conversion for 
unmodified and phosphorous modified CrTi samples. 
Temperature = 653 K; C 3 H 8 : O 2 = 2:1. 
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Propene Selectivity(%) 



Figure-4.11: Predicted propene selectivity versus propane conversion 
for unmodified and phosphorous modified CrTi samples. 
Temperature = 653 K; CsHg : O 2 = 2:1. 
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Chapter: 5 

Conclusions and Recommendations 

5.1. Conclusions 

The effect of phosphorus modification on the structure and reactivity of the 3% 
Cr203/Ti02 sample was studied by initially synthesizing four catalysts, by applying 
various characterization techniques (surface area, XRD, TPR and EPR) and finally 
carrying out the ODH of propane reaction. 

The surface area and XRD results revealed that the Ti02 support structure is not 
significantly affected with deposition of chromia and/or phosphorus. In the XRD analysis 
no new crystalline phase was detected suggesting that no new compounds were formed. 
The TPR studies gave two peaks due to the reduction of surface chromate group. The first 
reduction peak at ~615K corresponds to the reduction of the surface Cr^'*' species. The 
peak area associated with the surface chromium oxide species decreased with the 
phosphorus addition suggesting that the surface Cr^'^ species is progressively poisoned. 
Peaks due to y-Cr^'^ and clustered Cr^'^-O-Cr^'^ in amorphous CrP 04 were observed in the 
EPR spectra. The amount of Cr^"^ due to the amorphous CrP04 increases with phosphorus 
addition. 

The ODH of propane was performed with unmodified and phosphorus modified 
CrTi samples to examine its catalytic behavior. The propene selectivity was observed to 
increase and reach a maximum for a Cr:P ratio of 1.5 and then decrease with phosphorus 
modification. The conversion, however, decreased monotonically with phosphorus 
modification. 
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Two power law models were considered to explain the experimental data. The 
best suited power law model suggested propene as a primaiy product and carbon oxides 
as secondary product formed by propene degradation. The pre-exponential factors 
obtained for this model decreased with phosphorus addition (except for for ICrlPTi) 

at different rates, which was due to the decrease in number of active sites. The activation 
energies for the propene degradation to CO and CO 2 were found to be less than that of 
propane ODH. No specific trend in partial pressure exponents is observed as an effect of 
phosphorus addition to Cr 203 /Ti 02 catalyst. Based on the kinetic-parameters the propene 
jdeld was predicted and observed to be highest for l.SCrlPTi catalyst, which possessed 
the lowest (k 2 +k 3 )lki value. Thus, controlling the by changing the catalyst 

composition wall give rise to increased propene yields. 

5.2. Recommendations 

Based on the observations and conclusions of the present study, the following 
recommendations can be made; 

1. More detailed characterization of the modified catalysts using other techniques, such 
as XPS, DRS, dehydrated Raman spectroscopy and electron microscopy should be 
performed to obtain better insight into the nature of the catalysts. 

2. Preparation of the modified catalysts should be in various sequences, i. e., deposition 
of the active species to the modified-support or deposition of the modifier 
(phosphorus) to the supported catalysts, in order to understand the effect of 
preparation sequence on the structure-reactivity of the modified catalysts. 
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3. Mechanistic models, such as Langmuir-Knselwood and Eley-Rideal and Redox 
models should also be carried out for predicting the reaction pathways of propane to 
propene and for the better understanding of the reaction mechanism. 

4. Use of other modifiers to observe if the propene yield can be increased further. 
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Appendix-1 


Raw Peak Areas obtained in Contact Time Studies 

Table A- 1.1 Raw Peak Areas obtained in the contact time study of 3CrTi catalyst. 

Weight of the catalyst = 0. 1 g; Temperature = 653 K; Propane; Air=2:l . 


Total Flow 
Rate 
(cc/min) 

Contact 
Time (kg 
m^s) 

Raw Areas 

CsHg . 

C3H6 


CO 

120 

50 

1813.843 

2.989 



120 

50 

1815.875 

2.920 

4.269 

1.877 

90 

70 

1875.932 

_ 

3.542 

5.205 

2.309 

90 

70 

1874.602 

3.500 

5.387 

2.480 

75 

80 

1789.086 

3.761 

6.353 

2.491 

75 

80 

1756.879 



2.523 

45 

130 

1770.500 

4.255 

9.131 

3.354 

45 

130 

1774.844 

4.520 

8.302 

3.301 

30 


1637.594 

5.007 

10.955 

4.105 

30 

200 

1744.080 



4.981 

20 

300 

1697.442 

6.576 

15.250 

5.637 

20 

300 

1841.882 

7.078 

17.011 

6.027 
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Table A-1.2 Raw Peak Areas obtained in the contact time study of 2CrlPTi catalyst. 

Weight of the catalyst = 0.2 g; Temperature = 653 K; Propane; Air =2:1. 


Total 

Flow 

Rate 

(cc/min) 

Contact 
Time (kg 
m'^s) 

Raw Areas 

CaHs . 

CjHfi 

CO2 

CO 

75 

160 

1465.550 

3.471 

3.342 

1.412 

75 

■ 

160 

1460.524 

3.321 

2.989 

1.012 

60 

200 

1438.330 

4.085 

3.660 

1.740 

60 

200 


3.835 

3.577 

1.684 

45 

270 

1385.420 

4.561 

4.461 

1.942 

45 

270 

1350.102 

4.213 

4.587 

2.014 

30 

400 

1345.750 

5.675 

5.624 

2.831 

30 

400 

1330.425 

5.560 

5.501 

2.534 
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Table A-L3 Raw Peak Areas obtained in the contact time study of l.SCrlPTi catalyst. 

Weight of the catalyst = 0.2 g; Temperature = 653 K; Propane; Air =2:1. 


Total Flow 
Rate 
(cc/min) 

Contact 
Time (kg 
m'^s) 

Raw Areas 

C 3 H 8 

C 3 H 6 

CO 2 


120 

100 

1355.740 

5.157 

2.286 

1.371 

120 

100 

1423.050 

3.510 

1.730 

0.886 

90 

130 

1473.610 

3.990 

1.969 

1.039 

90 



4.010 

2.204 

1.222 

75 

160 

1515.892 

3.943 

1.986 

1.112 

75 

160 

1463.640 

3.676 

1.981 

1.127 

45 

270 

1463.790 




45 

270 

1448.164 

4.851 

2.847 

1.525 

30 

400 

1404.390 

5.719 

3.642 

1.874 

30 

400 

1418.528 

6.609 

4.189 

2.195 

20 

600 

1371.710 

8.531 

5.462 

2.683 

20 

600 

1371.836 

8.362 

5.208 

2.728 


80 





























































Table A-1.4 Raw Peak Areas obtained in the contact time study of ICrlPTi catalyst. 
Weight of the catalyst = 0.2 g; Temperature = 653 K; Propane: Air =2: 1 . 


Total Flow 
Rate 
(cc/min) 

Contact 
Time (kg 
m'^s) 

Raw Areas 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

120 

100 

1584.510 

2.481 

1.374 

0.924 

120 

100 

1560.964 

2.670 

2.170 

1.028 

90 

130 


2.910 

2.018 

1.158 

90 

130 

1540.741 

3.141 

1.654 

0.954 

75 

160 

1545.165 




75 



3.536 

2.453 

1.258 

45 

270 

1494.786 

4.148 

3.268 

1.643 

45 

270 

1500.742 

4.056 

3.158 

1.548 

30 

400 

1582.866 

5.165 

3.641 

1.741 

30 

400 

1522.866 

4.889 

3.841 

2.048 

20 

600 

1535.718 

5.904 

2.018 

1.158 

20 ■ 


1500.982 

5.748 

1.985 

1.425 
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Appendix-2 

Table A-2.1 Raw Peak Areas obtained in Propane ODH reaction study for 


3CrTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75cc/min. 


C 3 H 8/02 

ratio 

Zero 

Area 

Reaction 

Temperature 

(K) 

Raw Areas 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

1:1 

555 

633 




mm 

633 




1.843 

653 





653 





673 





673 





2:1 

1280 

633 

1255.146 

3,824 

22.904 

3.356 

633 

1260.799 

3.542 

20.394 

3.174 

653 

1254.294 




653 

1260.453 

9.338 

41.002 


673 

1260.831 



11.934 

673 

1273.408 

16.795 

48.213 

11.510 

3:1 

1442 

633 

1405.301 

7.317 

29.389 

7.025 

633 

1427.512 

4.940 

20.948 

6.681 

653 

1426.402 

27.087 

56.658 

19.368 

653 

1410.054 

23.226 

43.326 

16.633 

673 

1372.779 

34.176 

57.567 

23.607 

673 

1362.297 

40.128 

56.889 

26.231 


* The values used in the study are marked as bold. 
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Table A-2.2 Raw Peak Areas obtained in Propane ODH reaction study for 


2CrlPTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75cc/min. 


C 3 H 8/02 

ratio) 

Zero 

Area 

Reaction 

Temperature 

(K) 

Raw Areas 

C 3 H 8 

CsHfi 

CO 2 

CO 

1:1 

340 

633 

327.838 

0.785 

2.341 

0.649 

633 

338.846 

0.761 

2.179 

0.611 

653 




1.06 

653 

338.797 

1.453 

3.953 

1.077 

673 

331.327 

2.451 

5.729 

1.763 

673 




1.778 

2:1 

1300 

i 

633 

1280.757 

2.929 

5.652 

1.934 

633 





653 


7.685 



653 

1270.368 

7.801 

15.605 

4.623 

673 



31.797 

10.327 

673 

1236.435 

14.912 

23.816 


3:1 

— 

1450 

633 

1437.919 

3.195 

4.635 

1.708 

633 

1424.675 



1.682 

653 

1439.016 

6.398 

9.995 

3.542 

653 

1425.007 

6.600 

10.046 


673 

1415.524 

17.232 

25.148 

9.136 

673 

1413.159 

16.901 

23.028 

8.286 


* The values used in the study are marked as bold. 
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Table A-2.3 Raw Peak Areas obtained in Propane ODH reaction study for 


1 . 5Cr 1 PT i catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75cc/min. 


C 3 H 8/02 

ratio) 

Zero 

Area 

Reaction 

Temperature 

(K) 

Raw Areas 

CaHg 

C 3 H 6 

CO 2 

mm 


760 

633 

751.730 

1.333 

3.132 

mm 

633 

758.900 

1.204 

3.275 

0.957 

653 

741.316 

2.663 

5.967 

1.778 

653 

751.404 

2.730 

8.108 

1.955 


617.690 

3.924 

9.990 

2.994 

673 

610.255 

3.320 

8.698 


2;1 

1080 

633 

1045.534 

2.311 

3.433 

1.029 

633 

1044.583 



1.205 

653 




|[[^^| 

653 

IQH^j 

3.918 

5.848 

1.976 

673 

1022.888 

9.009 


IQQI 

673 




4.448 

3:1 

1960 

633 

1888.488 

4.529 

5.353 

2.088 

633 

1949.125 

4.926 

5.799 


653 

1846.969 

8.697 

11.333 

4.002 

653 




4.313 

673 




7.712 

673 

1823.903 

18.108 

23.367 

8.416 


* The values used in the study are marked as bold. 
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Table A-2.4 Raw Peak Areas obtained in Propane ODH reaction study for 
ICrlPTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75cc/min. 


C 3 H 8/02 

ratio) 

Zero 

Area 

Reaction 

Temperature 

(K) 

Raw Areas 

CaHg 

C 3 H 6 

CO 2 

mm 

1:1 

600 

633 

595.735 

1.542 

3.624 

0.384 

633 

588.081 

1.414 

3.525 

0.37 

653 

585.086 


6.762 

0.411 

653 

589.716 

2.290 

6.409 


673 

579.551 

2.966 


0.747 

673 

597.328 

2.867 


0.722 

, 

2:1 

765 

633 

751.151 



0.265 

633 

j 

748.965 




653 

734.224 




653 

759.245 




673 

726.647 

2.921 

2.118 

0.969 

673 

712.654 

2.883 



3:1 

1950 

633 

1945.708 

1.516 

1.389 

0.419 

633 

1925.895 

1.305 

1.312 

0.435 

653 

1880.649 

2.495 

2.149 

0.802 

653 

1881.576 

2.578 

1.924 

0.715 

673 

1732.371 

4.417 

3.173 

1.355 

673 

1711.762 

4.66 

3.467 

1.361 


* The values used in the study are marked as bold. 
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Appendix-3 

wmmmmKMmmmmmmmmmmmmmmmmmmrnmmmmmmmmmmm 


TABLE A-3.1: Input and output mole percentages for 3CrTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75 cc/min 


Temperature 

Input mole (%) 

Output mole (%) 

(K) 

C 3 H 8 


C 3 H 8 

C3H6 

CO 2 

CO 

633 

17.355 

17.355 

17.051 

0.049 

0.208 

0.042 

653 

17.355 

17.355 

16.229 

0.153 

0.504 

0.112 

673 

17.355 



17.355 

15.204 

0.331 

0.833 

0.283 

633 

29.577 

14.789 

29.289 

0.039 

0.322 

0.047 

653 

29.577 

14.789 

28.876 


0.597 

0.109 

673 

29.577 

14.789 

28.729 

0.228 

0.732 

0.192 

633 

38.650 

12.883 

38.363 

0.034 

0.231 

0.040 

653 

38.650 

12.883 

38.126 



0.081 

673 

38.650 

12.883 

37.062 



IIIQQH 
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TABLE-A-3.2: Input and output mole percentages for 2CrlPTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75 cc/min 


Temperature 

Input mole (%) 

Output mole (%) 

(K) 

C 3 H 8 

O 2 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

633 

17.355 

17.355 

17.724 

0.042 

0.104 

0.029 

653 

17.355 

17.355 

17.439 

0.078 

0.180 

0.051 

673 

17.355 

17.355 

17.330 

0.134 

0.274 

0.084 

633 




0.082 

0.121 

0.042 

653 

29.577 



0.185 

0.317 

0.097 

673 

29.577 

14.789 

28.349 

0.356 

0.500 

0.17 

633 


12.883 

38.470 

0.089 

0.114 

0.042 

653 

38.650 

12.883 

38.124 

0.185 

0.245 


673 

38.650 

12.883 

37.871 

0.483 

0.616 

0.224 
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TABLE-A-3.3: Input and output mole percentages for l.SCrlPTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75 cc/min 


Temperature 

Input mole (%) 

Output mole (%) 

(K) 

C 3 H 8 

O 2 

C 3 H 8 

C 3 H 6 

CO 2 

m 

633 

17.355 

17.355 


0.032 

0.066 

0.016 

653 

17.355 

17.355 




m 

673 

17.355 

17.355 





633 

29.577 

14.789 

29.013 

0.067 

0.087 

0.026 

653 

29.577 

14.789 

28.663 



0.084 

0.092 

0.039 

673 




0.262 



633 

38.650 

12.883 




0.039 

653 

38.650 

12.883 





673 

38.650 

12.883 



0.407 

0.144 
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TABLE-A-3.4: Input and output mole percentages for ICrlPTi catalyst 

Weight of the catalyst = 0.2 g; Total flow rate =75 cc/min 


Temperature 

Input mole (%) 

Output mole (%) 

(K) 

CaHs 

O 2 

C 3 H 8 

C 3 H 6 

CO 2 

CO 

633 

17.355 

17.355 

17.129 

0.047 

0.096 

0.009 

653 

17.355 

17.355 

17.085 

0.070 

0.170 

0.012 

673 

17.355 

17.355 

17.080 

0.090 

0.328 

0.020 

633 

29.577 

14.789 


0.113 

0.069 

0.032 

653 

29.577 

14.789 

27.818 

0.072 

0.054 

m 

673 

29.577 



0.114 

0.072 

0.033 

633 

38.650 

12.883 

38.248 

0.029 

0.022 


653 

38.650 


37.288 

0.052 

0.033 

0.011 

673 

38.650 

12.883 

35.521 

0.093 

0.059 

0.025 
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Appendix-4 


Table 4.7: Kinetic-parameters for PL-2 Model 

Calculated at 0-2958 atm, 0.148 atm, T„ = 653.16 K 


Parameter 

Units 

Catalyst 

3CrTi 

2CrlPTi 

l.SCrlPTi 

ICrlPTi 

^10 

“T 

s 

Co i-N 

S ^ 

0.99 

0.69 

0.50 

0.18 

^20 

0.33 

f 



0.03 

^30 

2.16 

0.47 

0.22 

0.26 

El 

-“H 

o 

s 

183 

140 

142 

110 

E2 


125 

117 

97 

Es 

106 



102 

124 

130 

ai 

Dimensionless 

0.28 

1.0 

1.0 


bi 

1.0 

0.35 

0.01 

1.0 

% 

0.01 

0.6 

1.0 

0.26 

bz 

0.27 

0.01 

1.0 

1.0 

% 

0.22 

0.65 

0.49 

0.01 

bs 

1.0 

1.0 

0.40 

1.0 
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